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The advent of solid state nanodevices allows for interrogating the physico-chemical properties of a 
polyelectrolyte chain by electrophoretically driving it through a nanopore. Salient dynamical aspects 
of the translocation process have been recently characterized by theoretical and computational 
studies of model polymer chains free from self-entanglement. However, sufficiently long equilibrated 
chains are necessarily knotted. The impact of such topological "defects" on the translocation process 
is largely unexplored, and is addressed in this study. By using Brownian dynamics simulations on a 
coarse-grained polyelectrolyte model we show that knots, despite being trapped at the pore entrance, 
do not per se cause the translocation process to jam. Rather, knots introduce an effective friction 
that increases with the applied force, and practically halts the translocation above a threshold force. 
The predicted dynamical crossover, which is experimentally verifiable, is of relevance in applicative 
contexts, such as DNA nanopore sequencing. 
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Nanopores, namely holes of nanoscale dimensions 
carved out of biological or solid-state membranes, are in- 
creasingly becoming an important tool to probe chemical 
and physical properties of polymers [THS]. For instance, 
polyelectrolytes, such as DNA, can be electrophoretically 
translocated through a pore and their chemical composi- 
tion can be inferred either by the ionic current blockaded 
by the DNA strand [7 , or from the electrical current 
measured perpendicular to the DNA backbone [8 . These 
approaches are being actively investigated because they 
hold great promise for fast and low-cost DNA sequencing. 

However, one of the issues that has received much less 
attention is related to possible limitations arising from 
the maximum length of the polymer that can be elec- 
trophoretically translocated through a nanopore without 
obstruction from the inevitable chain self-entanglement 
(knots). In fact, it is well known that, the incidence 
of knots increases exponentially with the chain contour 
length [9 , and, in turn, can affect kinetic, mechanical and 
equilibrium properties of sufficiently long (bio-) polymers 
P^QHTQ] . A knot in the polymer chain is then an unwanted 
potential obstruction to its translocation through a pore; 
much like the knot we customarily make at the end of a 
sewing thread to prevent it from "translocating" through 
a threaded piece of cloth. 

Motivated by these observations, in this Letter we set 
out to investigate theoretically and numerically the dy- 
namics of pore translocation of chains whose contour 
length exceeds by orders of magnitude their persistence 
length. On the one hand, this allows us to push to un- 
precedented contour lengths the assessment of the va- 
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FIG. 1: [Color online] (a) Snapshot of an instantaneous 
configuration of the 15/im-long model polyelectrolyte driven 
electrophoretically through the nanopore. The arrow indi- 
cates the translocation direction, (b) Cut-through view of the 
pore region for configuration in panel (a). The chain back- 
bone, with a tightened trefoil-knot at the pore entrance, is 
highlighted in the inset. 



lidity of previously-suggested dynamical scaling relation- 
ships for the translocation process of unknotted chains 
(see, e.g., Ref. pOHM] ). On the other hand, we can 
clarify the impact of spontaneous knotting on the driven 
translocation of biopolymers. This avenue appears to be 
largely unexplored except for the recent protein-related 
investigation of Huang and Makarov [18 . 

In particular, we show that knots do not per se cause 
the translocation process to halt. More precisely, they are 
found to act as plug-like obstructions of the pore only if a 
threshold driving force is exceeded. Based on this result 
it is expected that accounting for the topology-dependent 
dynamical crossover ought to be important for applica- 
tions that employ nanopores, such as the detection and 
sequencing of DNA filaments [U [2] . 
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We consider a model polyelectrolyte chain that is elec- 
trophoreticahy driven through a pore embedded in the 
slab separating the cis and trans semi-spaces [1], see Fig. 
[l] For definiteness, the salient physical properties of 
the (otherwise general) model system are set to match 
those of 15/im-long ssDNA filaments in a solution with 
O.IM monovalent salt and translocated through an arti- 
ficial nanopore. The nominal slab thickness and effective 
pore diameter were set respectively equal to lOnm and 
2nm, consistent with the typical geometry of solid-state 
nanopores (which are about 3-4 times longer than bio- 
logical nanopores) [1 . 

The polyelectrolyte chain consists of TV = 15000 beads, 
with diameter equal to the nominal ssDNA thickness, 
a = Inm. Because the monovalent counterions reduce 
the phosphates electrostatic charge by ~ 50% [27 the 
effective charge of the beads, each spanning three nu- 
cleotides, is set equal to q = — 1.5e. The chain potential 
energy is accordingly given by: 
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where dij is the distance of monomers i and j, and the 
three terms, which enforce respectively the chain connec- 
tivity constraint (/7fene), the pairwise Lennard- Jones in- 
teraction (/7lj) and Debye-Hueckel (/7dh) repulsion have 
the form: 
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where Rq = 1.5<j, /^b is the Boltzmann constant, T = 
300K, £ = 80 is the water dielectric constant, Adh = 
0.9nm is the Debye screening length, and Cij is equal to 
1 if |z — j| = 1, and 10 otherwise. 

A Monte Carlo scheme, employing unrestricted, non 
topology-preserving local and global moves [28l [29] , was 
first used to generate an equilibrated set of conforma- 
tions for flexible self- avoiding chains of N = 15000 beads 
in bulk. The degree and type of entanglement of the fil- 
aments was established using the minimally-interfering 
closure scheme [30l [31]. It was found that ~ 2% of the 
configurations were knotted and, in more than 90% of the 
cases, they consisted of the simplest knot type: the trefoil 
or 3i knot (see inset in Fig. for a trefoil knot repre- 
sentation). The knots spanned, on average, about 10% of 
the chain; the sizable knot length is consistent with the 
expected delocalisation of knots in unconstrained chains 



Since the main focus of the Letter concerns the impact 
of chain topology on the translocation dynamics, we ne- 
glect here the otherwise important issue of how a chain 
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FIG. 2: [Color online] (a) Average time, (r), required to 
translocate a fraction x of an unknotted chain at various 
driving forces, /. The average is taken over 20 uncorrelated 
unknotted conformations. Inset: statistical distribution of 
X at three different times and / = 20pN. Because the pro- 
cess timescale falls in the overdamped regime, time is ex- 
pressed in units of the nominal monomer self-diffusion time: 
Td — Tirja^ /2i^bT . For water viscosity, 77 = IcP, Td ^ 0.5ns 
[39^. (b) Collapse of the data-points in panel (a) using the 
theoretical dynamical scaling, see text. 



in the bulk approaches the pore and enters it [35]. Ac- 
cordingly, we extracted from the equilibrated ensemble 
several uncorrelated knotted and unknotted chains, for 
which a rigid global translation or rotation could bring 
one of the ends at the pore entrance while the chain re- 
mainder stays in the cis semispace. 

The beads inside the pore are driven through it by a 
constant force whose magnitude, /, typically chosen in 
the 4 — 40 pN range. Assuming that the electro-osmotic 
screening inside the channel j6l [36j reduces by ~ 50% 
the charge density of ssDNA, these forces correspond to 
a uniform electric field of 0.15 — 1.5V per lOnm acting on 
each bead. The translocation dynamics of the polyelec- 
trolyte chain is integrated numerically using the fixed- 
volume and constant-temperature molecular dynamics 
simulation scheme implemented in the LAMMPS package 
[37]. As in other coarse-grained approaches, no explicit 
hydrodynamic treatment is introduced [38] . 

Unknotted chain dynamics - The translocation dynam- 
ics for various pulling forces and averaged over an en- 
semble of 20 unknotted chains is illustrated in Fig. [2^. 
The dynamical process is sensitive to the specific geo- 
metrical arrangements of the pulled chains: in fact, the 
heterogeneity of the process increases with time. As 
shown in the inset at fixed force and at a given simu- 
lation time, the translocated fraction of the chain can 
range from 50% to 100% depending on the initial con- 
figuration. Notwithstanding these differences, the av- 
erage asymptotic translocation times appear to follow 
closely the dynamical scaling relationship predicted the- 
oretically, (r) (X x^^^ / f where x is the translocated frac- 
tion of the chain [20 -22 . In fact, by using the theoretical 
self- avoiding matrix exponent v = 0.59, we find the above 



3 



- •-• f=4pN 


1 mr 


f=10pN 




4^f=20pN 




^f=30pN 






1 1 



10000 1e+05 1e- 

time [xj 

FIG. 3: [Color online] Translocation kinetics of 
trefoil-knotted chain at various driving forces. 
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expression well satisfied, see Fig. 

It should also be noted that this scaling relationship, 
which was first formulated using dimensional and heuris- 
tic arguments [2QH22 . was more recently shown to hold 
only for asymptotically long chains [26^. In fact the non- 
equilibrium process that governs the propagation of the 
tensile disturbance along the chain is sensitive to finite 
chain effects [26^. The collapse of the data points in 
Fig. (2)3 clarify that the theoretically-predicted asymp- 
totic scaling relationship holds satisfactorily for chains of 
N = 15000 beads with the considered driving protocol. 

Knotted chain dynamics - Compared to the unknot- 
ted case, the translocation dynamics of knotted chains 
has instead a dramatic, non-monotonic dependence on 
the driving force. This is illustrated in Fig. [3] for one 
particular trefoil-knotted configuration where the knot 
occupies the mid-portion of the chain and spans about 
10% of the chain. It is seen that, at the smallest driv- 
ing force the chain translocation is well-consistent with 
the average translocation dynamics of unknotted chains. 
At higher forces, however, the standard dynamic scaling 
is satisfied only up to when about 50% of the chain is 
translocated, after which a noticeable slowing down of 
the process ensues. Notice that at the highest force, the 
translocation process appears to be practically halted. 

This dynamical crossover is understood by monitor- 
ing the instantaneous position and size of the knot along 
the chain. The data for / = 20pN are shown in Fig. |4] 
It is seen that the knot position along the chain is un- 
perturbed until it is reached by the propagating pulling 
front. At this stage, the progressive dragging of the chain 
through the pore causes the knot to tighten. Notice that, 
because the sequence index of the knot distal end remains 
about constant in time, the knot tightening process con- 
sists of the removal of the "slack" from the end that is 
nearest to the pore. After this stage, the tightened knot 
is localized at the pore entrance and the translocation 
proceeds by chain reptation through the knot "defect". 
These results illustrate the remarkable impact that non- 
trivial spontaneous entanglement of long polymer chains 
has on the translocation dynamics as a function of the 
pulling force. 



FIG. 4: [Color online] Dynamical evolution (a) of the translo- 
cated chain fraction and (b) of the relative knot location along 
the chain for the same translocation run at / = 20pN shown 
in Fig. [3] At time t ^ 7 lO^Td (dotted line) the propagating 
pulling front reaches the knot which progressively tightens 
until it reaches the pore at t ~ 1.2 IO^Tcz (dashed line). 



To quantify in the most transparent way the topology- 
dependent aspect of the effect, and separate it from the 
one associated to the chain geometry we have extended 
our analysis in two complementary directions. First, by 
comparing the dynamics of chains with same geometry, 
but different topology near the pore entrance. Secondly, 
by suitably averaging the translocation dynamics over 
chains with different geometry but same knot topology. 

For the first analysis, immediately after the knot is 
tightened at the pore entrance, one can locally perturb 
the chain geometry at the pore so as to untie the knot, 
while leaving unaltered the coordinates of all other parts 
of the chains, see Fig. [5^. Next, by following both the dy- 
namical evolution of the knotted and unknotted versions 
of the chain (with the same initial velocities for both sim- 
ulations) it is possible to compare the net effect of the 
localized knot defect on the chain translocation dynam- 
ics. The results are shown in Fig. [5] and aptly illustrate 
that the slowing down of the translocation process is as- 
cribable to the presence of the localized knot at the pore 
entrance. 

A quantitative assessment of the topology-dependent 
hindrance can be made by comparing the average translo- 
cation times of chains that are unknotted (see Fig. 2a) 
and of chains whose knot is initially close to the chain end 
at the pore entrance. Specifically, we considered chains 
where the knot was located within the first 10-20% of 
the chain, and for which the knot tightening and trap- 
ping occurs within a time span that is typically less than 
5% of the full translocation process. We computed the ef- 
fective friction coefficient, jknot^ of these knot-dominated 
processes and found that it has a dramatic dependence 
on /. This contrasts with the unentangled case where 70 
is practically force-independent, as testified by the col- 
lapse of the data in Fig. [2] The effect is illustrated in 
Fig-i 
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FIG. 5: [Color online] Dependence of the translocation dy- 
namics on the chain topology. After unknotting the chain by 
the shown local modification of the chain at the pore entrance 
the translocation velocity is dramatically enhanced. For clar- 
ity, here and in Fig. [6] the chain centerline is rendered as a 
continuous tube with diameter smaller than a. 
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FIG. 6: [Color online] Force dependence of the normalised in- 
verse friction coefficient, 7o/7fcnot for chains hosting 3i knots 
and 4i knots. At each force, 70 is defined as 70 = (rt)//cr, 
where (rt) is the average time required by unknotted chains 
to fully translocate. An analogous definition is used to com- 
pute 7fcnot for 3i- and 4 1 -knotted chains. Averages are taken 
over 20 configurations for the unknotted case and five config- 
urations for each non-trivial topology. 

Notice that at the lowest pulling force, / = 4pN, where 
a fairly tight (spanning 27 monomers on average) knot is 
trapped at the pore entrance, '^knot is statistically com- 
patible with 7o. We therefore conclude that knot local- 
ization at the pore entrance is not per se an impediment 



for translocation. It becomes so only when higher pulling 
forces cause the chain beads to interact very tightly, re- 
sulting in a rapidly increasing effective friction. Consis- 
tent with the dynamics in Figs. |4] and [5| the trend in 
Fig. |6] indicates that the translocation process is practi- 
cally halted at / - 30pN. 

It is interesting to notice that the observed impact of 
topology on dynamics differs from the case of knotted 
chains that are passively ejected out of a small spherical 
cavity [40l|41]. In such systems, knots - even when tight 
- reduce the ejection speed by only a factor of 2 — 3 [40j, 
possibly due to the relatively small magnitude of the force 
driving the spontaneous ejection. 

Yet, consistent with available numerical results for 
knot-controlled DNA ejection and protein transloca- 
tion dynamics [18 we do observe that at moderate driv- 
ing forces, the dynamical hindrance is appreciably higher 
for knots that are more complex (and rarer) than tre- 
foils. In fact, as it is shown in Fig. |6]the average inverse 
friction coefficient of five 4i -knotted chains is noticeably 
smaller than of 3i-knotted ones. Forces in the 30 — 40pN 
range are nevertheless sufficient to halt the translocation 
process of all such 4i-knotted configurations as well as of 
five instances of 5i- and 52-knotted ones (see supporting 
material) . 

Finally, we stress that the results presented here ought 
to be relevant in applicative contexts, such as genomic 
nanopore sequencing where the high-throughput demand 
pushes towards interrogating longer and longer unin- 
terrupted DNA filaments at pulling forces comparable 
to those considered in this work (TJ |2]. The high- 
throughput condition inevitably leads to significant chain 
self-entanglement, while high pulling forces may cause 
tight knots to halt the translocation process. It should 
also be pointed out that in actual ssDNA chains the latter 
effect will expectedly be more severe than in our model 
system because of both the ramified character of the 
molecule, and base-pairing effects. Since, the predicted 
force- dependent topological jamming is experimentally 
verifiable, we hope that the present investigation will mo- 
tivate further studies aimed at quantifying this effect. 
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SUPPORTING FIG. 1: Translocation kinetics at low and high force for polymer chains with three different knot topologies 
4i, 5i and 62 as shown in the insets. At low force, / = 4pN, the translocation proceeds steadily but it is stalled at / ~ 35pN 
for the 4i -knotted chain and / ~ 40pN for 5i- and 52-knotted chains. The jamming force shown in the examples is typical: 
5 instances of 4i knots were all jammed at / = 35pN and 5 instances of five-crossing knots were all jammed at / = 40pN. A 
time step as small as 0.003TLj,{rLj — cr^m/e is the Lennard- Jones time associated to monomer dynamics 39 ) was used to 
avoid chain breaking due to the build-up of tensile strain in jammed configurations. The default LAMMPS parametrization 
[39^ Tj = implies that the effective bead mass, m \s ^ 2.5 • 10~^^kg. We recall in fact that e = kbT {kb is the Boltzmann 
constant and T = 300K), a = Inm, and 77 = IcP. With this effective bead mass, tmd corresponds to ~ 6ns. By analysing the 
velocity- velocity autocorrelation function for internal chain beads, it is established that the overdamped Langevin regime sets 
in at time-scales ^ 0.5ns, i.e. comparable to Td- In our analysis we accordingly focused exclusively on the mass- independent 
MD evolution over time-spans longer than Td and have accordingly expressed time-scales in units of r^. 



